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ABSTRACT 

An optimization method is proposed to reduce the vibration of 
thin-plate structures. The method is based on a finite-element shell 
analysis, a modal analysis, and a structural optimization method. In 
the finite-element analysis, a triangular shell element with 18 degrees 
of freedom is used. In the optimization, the overall vibration energy 
of the structure is adopted as the objective function, and it is mini- 
mized at the given exciting frequency by varying the thickness of the 
elements. The technique of modal analysis is used to derive the sen- 
sitivity of the vibration energy with respect to the design variables. 
The sensitivity is represented by the sensitivities of both eigenvalues 
and eigenvectors. The optimum value is computed by the gradient 
projection method and a unidimensional search procedure under the 
constraint condition of constant weight. 

A computer code, based on the proposed method, is developed 
and is applied to design problems using a beam and a plate as test 
cases. It is confirmed that the vibration energy is reduced at the 
given exciting frequency. For the beam excited by a frequency 
slightly less than the fundamental natural frequency, the optimized 
shape is close to the beam of uniform strength. For the plate, the 
optimum shape is obtained such that the changes in thickness have 
the effect of adding a stiffener or a mass, 

INTRODUCTION 

The machine design process generally has several steps. A key 
step is the evaluation of a preliminary design to determine whether 
or not the planned performance features are achieved. A decision to 
change or approve the design is then made. Although these steps 
depend on the designer, an optimization technique can aid in obtain- 
ing a satisfactory solution. The process of design optimization con- 
sists of two stages: the first is the formulation of the design 
mathematically, for example, the minimization of weight. The sec- 
ond stage is the process of solving the problem. Considerable 
research has been done in these fields (Atrek, 1989), and the signi- 
ficance of optimum design has increased. In connection with analy- 
sis, the first author has reported the optimum statical design of a 
plate based on maximization of its bending rigidity (Inoue, 1988). 
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Takatsu (1991) and Lim et al. (1989) conducted research on the 
vibration of gear systems. The gearboxes used in these studies were 
made of steel plates. One technique that is often used to avoid 
resonance shifts the natural frequency for the purpose of controlling 
and reducing the vibration of structures like gearboxes. A mass or 
stiffener is added at some position on the structure to change the 
resonance. Although effective, it is difficult to select the position 
and to estimate a suitable shift of natural frequency. Another tech- 
nique to avoid resonance reduces the dynamic response, or transfer 
function. The dynamic response is used as the index of vibration in 
this technique. However, it is not always the best index for evalu- 
ating the overall vibration of the structure because it is usually 
evaluated at a specific point. The first author adopted the vibration 
energy as the index and proposed a method to design a plate (Inoue 
et al., 1990). In this research, the method was refined and a com- 
puter program was developed for the optimum design of thin-plate 
structures. The program primarily uses finite-element and modal 
analyses. The latter is used successfully to derive the sensitivity of 
the vibration energy with respect to the design variable. The pro- 
posed method is applied herein to some beam and plate design prob- 
lems in order to demonstrate the optimization process. The 
application to the design problem of shell structure is discussed in 
another paper (Inoue et a!., 1992). 

OPTIMIZATION BASED ON THE MINIMIZATION OF 
VIBRATION ENERGY 

Modal Analysis of Forced Vibration 

An elastic structure is divided into finite elements and discre- 
tized to N degrees of freedom. The equation of motion of the struc- 
ture is given in matrix form as 

[M]{u(t)} + [C]{u(t)} + [K]{u(t)} = W 


where |Mj is the mass matrix, [C| is the damping matrix, [K] is 
the stiffness matrix, {u(t)} is the displacement vector, and {f(t)} 
is the exciting force vector. When proportional damping is assumed, 
the damping matrix is represented by the linear combination of the ’ 
mass matrix and the stiffness matrix 


1 


[Cl = <*[M) + flK] (2) 


where a and /} are the coefficients that define the proportional 
damping matrix. The solutions of the eigenvalue problem are 
obtained easily from the eigenvalues and eigenvectors of the 
undamped vibration problem with the same matrices [M] and [K] 
The eigenvector has orthogonality in a broad sense. Let U\ 
be the eigenvector normalised with respect to mass matrix 

(*r} T [Mj{*,} = S„ (3) 

*jj ere ^rs * s Kronecker’s delta, and r and s define the r th and 
s mode. The modal matrix (*] is defined by the assembly of 
these eigenvectors 


1*1 = K*i>. .... {(6„>] (4) 


where n indicates the number of eigenvectors considered. From the 
normalized orthogonality of eigenvectors, any displacements (u(t)} 
can be represented approximately by their linear combination 


(u(t)} = £ f r (t){^ r } = l$]{f{t)} (5) 

r=l 


{«i(t)} = ju($j{f} e i"‘ ( 10 ) 

From Eqs. (9) and (10), the vector {f} is given by 

= (“> 

Therefore, the real and imaginary parts of r th component are 
obtained as follows: 

(« + ^r){^r} T {F} 

l(A r - <u 2 ) 2 + (« + PX T ) 2 

w 

~( A r ~ wJ ){^ r } T {F} 

OJ 

— ( A r - «*)* + (<* + /JA r ) 2 
w 

where R and I are the real and imaginary parts and {. . .} T 
denotes the transpose matrix. 

The vibration energy T of the structure is given by 


£rR - 

frl ~ 


When the structure concerned is excited by a harmonic force of 
angular frequency 

(f(t)} = {F}e^‘ (6) 


where {F} is the exiting force amplitude. The vibration of the 
structure is written by the same form 

MO) = = [*]{f}e Jwt (7) 


t = 2(0 ’} t [m]{o} 

= I{r} T I*l T [M][*l^ [M 

where ({/*} and {£*} show the conjugate complex of (C/J and 

{$)• The energy is, therefore, represented by the following 
expression: 


where {U} indicates the amplitude and is represented by the modal 

matrix. Substituting Eqs. (6) and (7) into Eq. (1) and premultl- , n , , 

plying it by ft] ,n independent algebraic equations are obtained T = _ £ (f B + 4 T | 

to evaluate the vector {£}. The r 6f “ component is given in the 2 '-1 ' 1 

following form by solving the equation: 


fr = 




" 2 ) + jw(« + £A r ) 


(r = l,2,...,n) 


( 8 ) 


In the expression, A r denotes the r th eigenvalue. Consequently, 
the amplitude {U} for the vibration displacement is obtained by 
substituting {{} into Eq. (7). 

Sensitivity of Vibration Energy 

The velocity (u(t)} can be expressed in a form similar to 
Eq. (5). The vibration velocity is given by 


m) = {0}e jwt = \^){(W ut 


(9) 


The sensitivity of the vibration energy with respect to the design 
variable Xj (i = l,...,m) is given by 


£T 


n 

E 

r=I 


dxz 


4i] 


(15) 


mce o, 0, w, and {F} in Eq. (12) are independent of the design 
variable, the sensitivity of the vibration energy is represented by the 
function of the sensitivity of the eigenvalue and the sensitivity of the 
eigenvector. The evaluation of the sensitivities of both eigenvalue 
and eigenvector is based on the method proposed by Fox and 
Kapoor (1968). 


However, the velocity can be derived directly by differentiating the 
displacement in Eq. (7): 
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Search for the Optimum Value 

The optimization design problem is expressed as follows: 

Minimize T({x}) 

Behavior constraint g({x}) = W — W({x}) = 0 

Side constraint ^({x}) = Xj - > 0 (16) 


The program was implemented on Digital Equipment Corpor- 
ation VAX computers at the NASA Lewis Research Center. It is 
executed interactively for small models. In the case of a larger 
model like a three-dimensional structure, the program is submitted 
to the Cray Research Corporation XMP 4/28 computer. The pro- 
gram is written in FORTRAN 77. 


APPLICATION TO THE OPTIMUM DESIGN OF BEAM AND 
h 2 ({x}) = - Xi > 0 PLATE 


Design variables {x} = (x l ,x 2 ,...,xJ T 


where W and W({x}) are the prescribed weight and the function 
of the structure weight, respectively. The behavior constraint gives 
the design of constant weight. The element thickness is adopted as 
the design variable. The side constraints limit the thickness within 
the given minimum and maximum. 

In the optimization process, the design variables are modified 
along the feasible direction vector {d^} to reduce the objective 
function, and the new design variables are obtained in the following 
form: 

{x (k+1) } = (x< k >} + s{d (k) } (17) 


The superscript (k) means the k th iteration in the optimization 
process, and s is the step. If the variables are in the feasible 
region, the feasible direction vector is given by the steepest descent 
vector and is composed of the sensitivity of vibration energy with a 
minus sign. After the design variables reach the intersection formed 
by the constraint surfaces, the vector is projected on the inter- 
section. This procedure is the gradient projection method (Rosen, 
1961). The method can be applied to problems with linear or non- 
linear constraints. Since the constraints for the design problem dis- 
cussed in this research are linear, the feasible direction is easily 
determined. In order to obtain the step s in Eq. (17), the region 
which includes the relative minimum of the energy is searched by an 
iterative interpolation process. Once the region is found, the ener- 
gies are evaluated at the limits and the middle point of the region. 
Next a quadratic expression is obtained so that it goes through these 
points. The step size is determined to obtain the minimum of the 
expression. The region cannot be found every time because the 
objective function frequently decreases monotonously. In this case, 
the search procedure is stopped if the minimum is not obtained after 
10 iterations. Then the step at the last iteration is used, the sensi- 
tivity is calculated again, and the optimum is searched. The value 
of the energy at the last iteration is stored and used to determine 
convergence. The convergence criterion is 1/1000 relative variation. 

OPTIMUM DESIGN PROGRAM 

A computer program was developed based on the procedures 
described in the previous sections. The program consists of two 
main parts. One part is the analyzer, which is composed of the 
finite-element method, the eigenvalue analysis, and the evaluation of 
the vibration energy by means of modal analysis. The other part is 
the optimizer which includes the sensitivity analysis and search pro- 
cedure for the optimum value. 

Several finite-element programs are available; however, in this 
research, the program is written to make the program as compact 
and flexible as possible. A triangular shell element with 18 degrees 
of freedom is used. The element is formed by the combination of 
elements for plate bending and plane problem analyses, and both are 
given in Zienkiewicz and Cheung (1967). The consistent mass 
matrix is evaluated by using the method of numerical integration. 


In order to demonstrate the optimization process, the proposed 
method is applied to the optimum design problem of beams and 
plates. These constants were used in this research: modulus of 
elasticity in tension ~ 206 GPa, Poisson’s ratio = 0.3, mass dens- 
ity “ 8000 kg/m 3 , and the coefficients for composing the propor- 
tional damping matrix are a = 1.0 sec -1 and ji = 5.0xl0~ 7 sec. 

Optimum Shape of Simply Supported Beams 

The element used in this research is not a beam element but a 
shell element, and it is not best suited for the analysis of beams. 
However, a sufficiently narrow-width plate may be approximated as 
a beam. The pattern of the mesh is shown in Fig. 1, The element 
numbers are from 1 to 40 as shown in the figure. The beam has 
dimensions of length L = 200 mm, width B = 10 mm, and the ini- 
tial thickness t Q — 5 mm. 

The calculated lower four natural frequencies and their errors 
(compared to the theoretical values for a simply supported beam) 
are 288.9 Hz (0.5 percent), 1172 Hz (1.9 percent), 2704 Hz 
(4.5 percent), and 4982 Hz (8.3 percent), respectively. When the 
beam center is excited by a unit harmonic force, the frequency 
response of the vibration energy is evaluated at every 10 Hz, based 
on these lower four modes, and is illustrated by the dashed curve in 
Fig. 2. The abscissa indicates the exciting frequency. The peak in 
the figure approximately represents the resonance at the fundamen- 
tal natural frequency of 288.9 Hz. If this beam is excited by a 
unit harmonic force (1 N) at 270 Hz, the energy level of about 
0.207 xlO” 3 J is considerably high because of the influence of the 
resonance. The optimum design program was executed to minimize 
the vibration energy at the resonant frequency. The upper and 
lower limits of the thickness are given as 1 < t < 10. The process of 
energy convergence, which is represented by the ratio of the initial 
energy, is shown in Fig. 3. The energy decreases rapidly and con- 
verges to about one-tenth of the initial energy after four iterations. 
The frequency response of the optimally designed beam is indicated 
by the solid curve in Fig. 2. The fundamental natural frequency 
shifts to 319.7 Hz. 

The optimum shape of the beam is illustrated in Fig. 4. The 
thickness is maximum at the center and decreases gradually toward 
the supported ends. The shape is compared with the shape of a 
beam of uniform strength in Fig. 5. The uniform-strength beam is 
well known as the optimum shape of the static maximum stiffness 
design (Huang, 1968). A similar design performed by using a fine 
mesh of 20 by 2 by 2 gave a shape approximately the same as the 
aforementioned optimum shape in Fig. 5. 

In the optimization, the thicknesses of four elements in every 
section are forced to be equal. If the constraint is eliminated, the 
optimum shape illustrated in Fig. 6 is obtained. Most of the ele- 
ments reach the higher and lower thickness limits. However, the 
tendency for the element thickness to decrease toward the support- 
ing ends is the same as in the previous example. The element thick- 
ness variation is summarized in Table 1. The element thickness in 
each section is very close to each other in the first few steps of the 
process and then the difference increases gradually. The frequency 
response and the energy convergence are shown in Figs. 7 and 8, 
respectively. Compared with the former example, more iterations 
are needed until the energy converges, and the energy greatly 
reduces. The beneficial effect of this method on the reduction of 
vibration is confirmed by the vibration amplitude shown in Fig. 9. 
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Figure 1 . — Pattern of mesh used for the optimum design of 
beam. Length L = 200 mm, width B = 10 mm, initial 
thickness t Q = 5 mm. 



gure 2.— Frequency response of simply supported beam due 
to unit exciting force. Excitation frequency f ex = 270 Hz. 



Figure 3.— Convergence of vibration energy. 


Figure 4. — Optimum shape of simply supported beam. Excitation 
frequency f ex = 270 Hz. (The thicknesses of four elements in 
every section are forced to be equal.) 



Figure 5.— Comparison of the optimum shape with the shape of 
the beam of uniform strength. 



Figure 6. — Optimum shape of simp/y supported beam. Excitation 
frequency f ex = 270 Hz. 



Energy, J 


TABLE 1.— VARIATION OF ELEMENT THICKNESS OF SIMPLY SUPPORTED BEAM 



Frequency, Hz Number of iterations 

Figure 7. Frequency response (due to unit exciting force) of Figure 8. — Convergence of vibration energy with variable thick 

simply supported beam with variable thickness and width. ness and width. 

Excitation frequency, f ex = 270 Hz. 




.6x1 CH 



Figure 9.— Displacements due to continuative vibration (initials; 
uniform thickness, optimum 1, beam in fig. 4, optimum 2, 
beam in fig. 6). 


Optimum Shape of Simply Supported Plate 

The plate (whose dimensions are a = 600 mm, b =s 400 mm, 
and initial thickness t 0 = 5 mm) is divided into 6 by 2 by 2 
elements as shown in Fig. 10. The thickness constraint is given as 
2 5s t < 10. The eigenvalue problem Is solved first, and the lower 
12 modes are used for the optimization because they fully cover the 
frequency response range considered in the optimization. 

In the first example, the plate is excited by a central unit 
harmonic force (1 N) with a frequency of 370 Hz, which is slightly 
lower than the natural frequency of mode (3,1), 376 Hz. Mode (3,1) 
means the mode shape has two nodal lines parallel to the y-axis. 

The frequency response of the plate with initial thickness is eval- 
uated at every 10 Hz up to 1000 Hz and is shown by the dashed 
curve in Fi|. 11. The energy at the exciting frequency is about 
0.338 x10 J. The energy is considerably high because of the reson- 
ance influence. The energy convergence process is shown in Fig. 12. 
The energy decreases sharply by the first optimization and converges 
to about one- thousandth of the initial energy after eight iterations. 
The final frequency response is also shown by the solid curve in 
Fig. 11. The resonance peak under consideration shifts to about 
520 Hz, which causes a large vibration reduction at the excitation 
frequency. The optimum shape of the plate is illustrated In Fig. 13. 
A marked change of the shape is the increase in thickness of the ele- 
ments along the longer centerline which may function as a stiffener 
and lead to the vibration reduction of mode (3,1). 

In the second example, the same plate is excited by 390 Hz, 
which is a little higher than the natural frequency under consider- 
ation. The solution converges to the optimum shape as illustrated 
in Fig. 14. The frequency response and the energy convergence are 
shown in Figs. 15 and 16, respectively. In this example, the thick- 
ness increases at the central portion of the plate and appears to 
function as a mass. Indeed, the natural frequency under considera- 
tion is reduced by the optimization and shifts to about 314 Hz. 



Figure 10. Pattern of mesh used for the optimum design of plate 
(a = 600 mm, b = 400 mm, initial thickness t 0 «5 mm). 



Frequency, Hz 

Figure 11. — Frequency response of simply supported plate due 
to unit exciting force. Excitation frequency f ex = 370 Hz. 


In comparison with the first example, the energy at 390 Hz is 
higher than the former energy at the same frequency, thus demon- 
strating that the minimum energy obtained in these examples is a 
local minimum. The comparison also indicates that the natural fre- 
quency shift to a higher frequency region is better than the shift to a 
lower frequency region for the purpose of reducing the vibration. 


6 


Ratio of vibration energy 


1 



Figure 12.— Convergence of vibration energy. 



Frequency, Hz 

Figure 1 5. — Frequency response of simply supported plate due 
to unit exciting force. Excitation frequency f ex = 390 Hz. 



Figure 14. — Optimum shape of simply supported plate. 
Excitation frequency f ex = 390 Hz. 


Figure 1 6. — Convergence of vibration energy. 


CONCLUSION 


The vibration energy was selected as the index for vibration 
evaluation for the purpose of designing low-vibration thin-plate 
structures. A method was proposed to minimize the energy. The 
method is based on a finite-element shell analysis, a modal analysis, 
and a structural optimization technique. The vibration energy sensi- 
tivity with respect to the design variable was derived and shown to 
be expressed by the sensitivities of both eigenvalues and 
eigenvectors. 

The computer program was developed and applied to the 
analysis of a beam and a plate to study the optimization process. 
This method was demonstrated to be effective in reducing the energy 
at the given frequency. In the design of a plate, the optimum shape 
was obtained such that the increase in thickness functioned as a stif- 
fener or a mass. If the optimization frequency is slightly lower than 
the natural frequency, then this method shifts the natural frequency 
under consideration to a higher frequency region. The opposite is 
also true. Consequently, this method can control the direction of 
shift of the natural frequency. 
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